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CN ! Abstract 
> 

. Q ' In wireless networks with random node distribution, the underlying point process model and the 

I 

channel fading process are usually considered separately. A unified framework is introduced that permits 

^ '. 

the geometric characterization of fading by incorporating the fading process into the point process model. 
Concretely, assuming nodes are distributed in a stationary Poisson point process in R'^, the properties 
of the point processes that describe the path loss with fading are analyzed. The main applications are 
connectivity and broadcasting. 
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I. Introduction and System Model 

A. Motivation 



The path loss over a wireless link is well modeled by the product of a distance component (often called 
large-scale path loss) and a fading component (called small-scale fading or shadowing). It is usually 
assumed that the distance part is deterministic while the fading part is modeled as a random process. 
This distinction, however, does not apply to many types of wireless networks, where the distance itself 
is subject to uncertainty. In this case it may be beneficial to consider the distance and fading uncertainty 
jointly, i.e., to define a stochastic point process that incorporates both. Equivalently, one may regard the 
distance uncertainty as a large-scale fading component and the multipath fading uncertainty as small-scale 
fading component. 

This paper is an extension of preliminary work that has appeared at ISIT 2006, Seattle, WA, and ISIT 2007, Nice, France. 
M. Haenggi is with the Department of Electrical Engineering, University of Notre Dame, Notre Dame, IN 46556, USA. E-mail: 

mhaenggi@nd . edu 



February 2, 2008 



DRAFT 



2 



We introduce a framework that offers such a geometrical interpretation of fading and some new insight 
into its effect on the network. To obtain concrete analytical results, we will often use the Nakagami-m 
fading model, which is fairly general and offers the advantage of including the special cases of Rayleigh 
fading and no fading for m = 1 and m ^ oo, respectively. 

The two main applications of the theoretical foundations laid in Section 2 are connectivity (Section 3) 
and broadcasting (Section 4). 

Connectivity. We characterize the geometric properties of the set of nodes that are directly connected 
to the origin for arbitrary fading models, generahzing the results in [1], [2]. We also show that if the 
path loss exponent equals the number of network dimension, any fading model (with unit mean) is 
distribution-preserving in a sense made precise later. 

Broadcasting. We are interested in the single-hop broadcast transport capacity, i.e., the cumulated 
distance-weighted rate summed over the set of nodes that can successfully decode a message sent from a 
transmitter at the origin. In particular, we prove that if the path loss exponent is smaller than the number 
of network dimensions plus one, this transport capacity can be made arbitrarily large by letting the rate 
of transmission approach 0. 

In Section 5, we discuss several other apphcations, including the maximum transmission distance, 
probabilistic progress, the effect of retransmissions, and localization. 

B. Notation and symbols 

For convenient reference, we provide a list of the symbols and variables used in the paper. Most of 
them are also explained in the text. Note that slanted sans-serif symbols such as x and f denote random 
variables, in contrast to x and / that are standard real numbers or "dummy" variables. Since we model 
the distribution of the network nodes as a stochastic point process, we use the terms points and nodes 
interchangeably. 
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Symbol 


Definition/explanation 


[k] 


the set {1,2,... ,k} 


1a{x) 


indicator function 


u{x) 


— ^{x^o}{x) (unit step function) 


d 


number of dimensions of the network 





origin in M'^ 


B 


a Borel subset of M or M*^ 


Cd 


= 7r'^/Vr(l + d/2) 




(volume of the d-dim. unit ball) 


a 


path loss exponent 


5 


= d/a 


A 


= (d + l)/a 


s 


minimum path gain for connection 


F,f 


fading distribution (cdf), fading rv. 


Fx 


distribution of random variable X (cdf) 


$ = {Xi} 


path loss process before fading (PLP) 


s = m 


path loss process with fading (PLPF) 


^ = {Xi} 


points in $ connected to origin 


s = m 


points in S connected to origin 


A, A 


counting measure and density for $ 


N = S(M+) 


number of nodes connected to o 


#A 


cardinaUty of A 



C. Poisson point process model 

A well accepted model for the node distribution in wireless networks' is the homogeneous Poisson 
point process (PPP) of intensity A. Without loss of generality, we can assume A = 1 (scale-invariance). 
Node distribution. Let the set {/j}, i G N consist of the points of a stationary Poisson point process in 
of intensity 1, ordered according to their EucUdean distance ||yi — o|| to the origin o. Define a new 

'in particular, if nodes move around randomly and independently, or if sensor nodes are deployed from an airplane in large 
quantities. 
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one-dimensional (generally inhomogeneous) PPP {/i — ||yi — o||} such that < ri < r2 < . . . a.s. Let 
a > be the path loss exponent of the network and $ = {xj = r"} be the path loss process (before 
fading) (PLP). Let {f ,fi,f2, . . .} be an iid stochastic process with f drawn from a distribution F = Ff 
with unit mean, i.e., ¥,f = 1, and supp c Finally, let E = {^i = Xi/fi} be the path loss process 
with fading (PLPF). In order to treat the case of no fading in the same framework, we will allow the 
degenerate case F{x) = u{x — 1), resulting in <1> = H. Note that the fading is static (unless mentioned 
otherwise), and that {^j} is no longer ordered in general. We will also interpret these point processes as 
random counting measures, e.g., $(-B) = #{$ PI B} for any Borel subset B of M. 

Connectivity. We are interested in connectivity to the origin. A node i is connected if its path loss is 
smaller than i.e., if < The processes of connected nodes are denoted as $ = {xj : < 
(PLP) and 3 = {ii: ii< 1/4 = S n [0, l/s) (PLPF). 

Counting measures. Let A be the counting measure associated with i.e., A{B) = K^{B) for Borel 
B. For A([0, a)) = E<I>([0,a)), we will also use the shortcut A(a). Similarly, let A be the counting 
measure for ^. All the point processes considered admit a density. Let X{x) = dA{x)/dx and and 
X{x) = dA(x)/dx be the densities of <I> and ^, respectively. 

Fading model. To obtain concrete results, we frequently use the Nakagami-m (power) fading model. 
The distribution and density are 

T{m) 

. , m"*a;"*~^ cxp(— ma;) 
1 (m) 

where Fjc denotes the upper incomplete gamma function. This distribution is a single -parameter version 
of the gamma distribution where both parameters are the same such that the mean is 1 always. 

D. The standard network 

For ease of exposition, we often consider a standard network^ that has the following parameters: 
5 = d/a = \ (path loss exponent equals the number of dimensions) and Rayleigh fading, i.e., F{x) = 
(1 — e~'-^)u{x). 

Fig. 1 shows a PPP of intensity 1 in a 16 x 16 square, with the nodes marked that can be reached 
from the center, assuming a path gain threshold of s = 0.1. The disk shows the maximum transmission 
distance in the non-fading case. 

^The term "standard" here refers to the fact that in this case the analytical expressions are particularly simple. We do not 
claim that these parameters are the ones most frequently observed in reality. 
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Fig. 1. A Poisson point process of intensity 1 in a 16 x 16 square. The reachable nodes by the center node are indicated by a 
bold X for a path gain threshold of s = 0.1, a path loss exponent of a = 2, and Rayleigh fading (standard network). The circle 
indicates the range of successful transmission in the non-fading case. Its radius is 1/^/s w 3.16, and there are about tt/s w 31 
nodes inside. 



II. Properties of the Point Processes 
Proposition 1 The processes S, and S are Poisson. 

Proof: {ui} is Poisson by definition, so {rj} and $ = {xi} are Poisson by the mapping theorem 
[3]. S is Poisson since ^ is iid, and H(M) = H([0, 1/s)). ■ 
The Poisson property of ^ will be established in Prop. 6. 

Cor. 2 states some basic facts about these point processes that result from their Poisson property. 

Corollary 2 (Basic properties.) 

(a) A(x) = E<I>([0, x)) = CdX^ and \{x) = Cd5x^~^. In particular, for 5 = 1,$ is stationary (on M.'^). 

(b) ri is governed by the generalized gamma pdf 
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and Xj is distributed according to the cdf 



^x.(x) = l-^^^^^^^,. (4) 



The expected path loss without fading is 



E^i = cr^-;T7^- (5) 



T{i) 

In particular, for the standard network, the Xj are Erlang with Exj = i/cd. 
(c) The distribution function of is 



For 6 = 1 and Nakagami-m fading, the pdf of (,i is 



In particular, 



and 



(8) 



Till 

Mi = form>\ (9) 

Cdim - 1) 

-r,- ^ m^iim + i — l) ^ ^ 

Var = 2r ^^2^ 1^ for m > 2 . (10) 



^^.(^) = (r^) • (11) 



For the standard networks, 

F.(x) = ( 

CdX + l 

Proof: 

(a) Since the original d-dimensional process {/j} is stationary, the expected number of points in a ball 
of radius x around the origin is Cdx'^. The one-dimensional process {ri\ has the same number of 
points in [0,aj), and Xj = r", so E$([0,x)) = CdX^. For 5 = 1, A(x) = Cd is constant. 

(b) Follows directly from the fact that {yi} is stationary Poisson. ((3) has been established in [4].) 

(c) The cdf P[^i < x] is 1 — Exj(-F'(xi/x)) with Xj distributed according to (4). (7) is obtained by 
straightforward (but tedious) calculation. 

■ 

Remarks: 

- For general (rational) values of m, d, and a, F^. can be expressed using hypergeometric functions. 
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- (8) approaches 1 — exp(— c^x) asm ^ oo, which is the distribution of xi. Similariy, limj„_>oo = 
i/cd = Exj and lim^^oo Var = i/c^ = Varxj. 

- Alternatively we could consider the path gain process Since F^-i{x) = 1 — F^.(l/x), the 
distribution functions look similar. 

- In the standard network, the expected path loss E^j does not exist for any i, and for i = I, the 
expected path gain is infinite, too, since both xi and f are exponentially distributed. For i > 1, 

= Cd/{i - 1), and for i > 2, Var(er') = 2c2/((i - - 2)). 

- For the standard network, the differential entropy h{^i) = E[— In (^j)] is 2 — log for i = 1 and 
grows logarithmically with i For Nakagami-m fading h{^i) = 1 + 1/m — log c^. For the path gain 
process in the standard network, the entropy has the simple expression 

/^(C^) = ^+log(-), (12) 
which is monotonically decreasing, reflecting the fact that the variance Var^~^ is decreasing with 

- The are not independent since the Xj are ordered. For example, in the case of the standard network, 
the difference x^+i — Xj is exponentially distributed with mean I/q, thus the joint pdf is 

/xi...x„(a;i, . . . ,x„) = cSe-'^<*^"lo<a:i<...<cc„ , (13) 
where lo<a;i<...<a;„ denotes the (positive) order cone (or hyperoctant) in n dimensions. 

Proposition 3 For 5 = 1 and any fading distribution F with mean 1, 

S(S) = $(B) VB c R+ , 

i.e., fading is distribution-preserving. 

Proof: Since S is Poisson, independence of S(Si) and 2(^2) for Si PI B2 = is guaranteed. So 
it remains to be shown that the intensities (or, equivalently, the counting measures on Borel sets) are the 
same. This is the case if for all a > 0, 

E (#{xi : Xi > a, < a}) = E (#{xi : Xj < a, .^i > «}) , 

i.e., the expected numbers of nodes crossing a from the left (leaving the interval [0, a)) and the right 
(entering the same interval) are equal. This condition can be expressed as 

i-a poo 

/ X{x)F{x/a)dx= / X{x){l-F{x/a))dx Va>0. 

Jo Ja 
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If (5 



= 1, X{x) = Cd, and the condition reduces to 




/ F{x)dx= / (1-F(x))d; 



IX 



Jo Jl 



which holds since 




Jo Jl 
^ ' 



An immediate consequence is that a receiver cannot decide on the amount of fading present in the 
network if 5 = 1 and geograpiiical distances are not known. 

Corollary 4 For Nakagami-m fading, 6 = 1, and any a > 0, the expected number of nodes with Xi < a 
and > a, i.e., nodes that leave the interval [0, a) due to fading, is 



The same number of nodes is expected to enter this interval. For Rayleigh fading (m = 1), the fraction 
of nodes leaving any interval [0, a) is 1/e. 

Proof: E (#{xi : Xj < a}) = A(a) = c^a, and for Nakagami-m, the fraction of nodes leaving the 
interval is 



Clearly, fading can be interpreted as a stochastic mapping from Xj to ^j. So, {xj} are the points in 
the geographical domain (they indicate distance), whereas {^j} are the points in the path loss domain, 
since ^ is the actual path loss including fading. This mapping results in a partial reordering of the 



Fig. 3 illustrates the situation for 200 nodes randomly chosen from [0, 5] with a threshold s = 1. Before 
fading, we expect 40 nodes inside. From these, a fraction is moving out (right triangles), the rest 
stays in (marked by x). From the ones outside, a fraction (1 — e~^){ae) « 9% moves in (left triangles), 
the rest stays out (circles). 

For the standard network, the probability of point reordering due to fading can be calculated explicitly. 



E (#{xj : Xj < a, > a}) = c^a 




—m 



(14) 





Let Pi J = > Ci+j]. By this definition. 




(15) 
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Fig. 2. The points of a Poisson point process Xi are mapped and reordered according to ~ Xi/fi, where fi is iid exponential 
with unit mean. In the lower axis, the nodes to the left of the threshold 1/s are connected to the origin (path loss smaller than 
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Fig. 3. Illustration of the Rayleigh mapping. 200 points Xi are chosen uniformly randomly in [0,5]. Plotted are the points 
{xi, Xi/fi), where the /, are drawn iid exponential with mean 1. Consider the interval [0, 1] (i.e., assume a threshold s = 1). 

Points marked by x are points that remain inside [0, 1], those marked by o remain outside, the ones marked with left- and 
right-pointing triangles are the ones that moved in and out, respectively. The node marked with a double triangle is the furthest 
reachable node. On average the same number of nodes move in and out. Note that not all points are shown, since a fraction 
is mapped outside of [0, 5]. 
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Xj is Erlang with parameters i and Cd, yj is the distance from Xj to Xj+j and thus Erlang with parameters 
j and Cd, and the cdf of z := fn/fn+m is Fz{x) = x/{x + 1). Hence 



Pi J does not depend on c^. Closed-form expressions include Pi,i = 1 — In 2 0.307, and Pi,2 = 
3 — 41n2 Si 0.227. Generally ^ can be determined analytically. For k = 1,2,3,4, we obtain 1 — 
In 2, 12 In 2 -8, 167/2-1201n2, 1120 In 2- 776. Further, linifc^oo Pfc,fc = 1/3, which is the probability 
that an exponential random variable is larger than another one that has twice the mean. 

In the hmit, as i ^ oo, Pij = + 1), which is the probabihty that a node has the largest fading 
coefficient among j + 1 nodes that are at the same distance. Indeed, as i ^ oo, Xj+j < Xj(l + e) a.s. for 
any e > and finite j. 

While the are dependent, it is often useful to consider a set of independent random variables, obtained 
by conditioning the process on having a certain number of nodes n in an interval [0, a) (or, equivalently, 
conditioning on Xn+i = a) and randomly permuting the n nodes. In doing so, the n points {xj} and 
i = 1, 2, . . . , n are iid distributed as follows. 

Corollary 5 Conditioned on x„_|_i = a: 

(a) The nodes {xj}"^j^ are iid distributed with 





^ X < a 



(16) 



and cdf F^^(x) = [x/aY- 
(b) The path loss with fading 



n 



is distributed as 




(17) 



(c) For the standard network, 




(18) 



(d) For Rayleigh fading and S 



1/2. 




(19) 



Proof: As in (6), the cdf is given by 1 — E{F{y/x)) with y distributed as (16). 
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III. Connectivity 
Here we investigate the processes $ and H = H n [0, of connected nodes. 

A. Single-transmission connectivity and fading gain 

Proposition 6 (Connectivity) Let a transmitter situated at the origin transmit a single message, and 
assume that nodes with path loss smaller than 1/s can decode, i.e., are connected. We have: 

(a) $ is Poisson with X{x) = A(x)(l — F{sx)). 

(b) With Nakagami-m fading, the number N = $(R"'") of connected nodes is Poisson with mean 



{■ms)° r(m) 



and the connectivity fading gain, defined as the ratio of the expected numbers of connected nodes 
with and without fading, is 



EiVoo rn^ r(m) 
Proof: 

(a) The effect of fading on the connectivity is independent (non-homogeneous) thinning by 1 — F{sx 
¥[x/f < 1/s]. 

(b) Using (a), the expected number of connected nodes is 



/ A(x)da:; = / CdSx° — ^— - — -dx 
Jo Jo r(m) 



which equals KNm in the assertion. Without fading, EiVoo = lim^^oo = M^/s) = CdS ^, wiiich 
results in the ratio (21). 

■ 

Remarks: 

1) (20) is a generalization of a result in [1] where the connectivity of a node in a two-dimensional 
network with Rayleigh fading was studied. 

2) EiV can also be expressed as 

oo 

EiV = ^P[ei<l/s]. (22) 

i=l 

The relationship with part (b) can be viewed as a simple instance of Campbell's theorem [5]. Since 
N is Poisson, the probability of isolation is ¥{N = 0) = exp(— E7V). 

3) ETVi = CdS~^V{5 + 1), and EiVoo = c^s"''. For 5 = 1, N does not depend on the type (or presence) 
of fading. 
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Fig. 4. Connectivity fading gain for Nakagami-m fading as a function of 5 G [0, 3/2] and m £ [1, 5]. For S — 1, the gain is 
1 independent of m (thick line). 



4) The connectivity fading gain equals the 6-th moment of the fading distribution, which, by definition, 
approaches one as the fading vanishes, i.e., as m ^ oo. For a fixed 6, it is decreasing in m if 6 > 1, 
increasing if 6 < I, and equal to 1 for all m if 6 = 1. It also equals 1 if 5 = 0. For a fixed m, 
it is not monotonic with 6, but exhibits a minimum at some 6„im G i^A)- The fading gain as a 
function of 6 and m is plotted in Fig. 4. For Rayleigh fading and 5 = 1/2, the fading gain is it/ 2, 
and the minimum is assumed at 6mm ~ 0.462, corresponding to a 4.33 for d = 2. So, depending 
on the type of fading and the ratio of the number of network dimensions to the path loss exponent 
a, fading can increase or decrease the number of connected nodes. 

5) For the standard network, EiV = q/s and the probability of isolation is e~^'^l^ . 

6) The expected number of connected nodes N°- with x, < a is 

^N" = Cda^Fl{\ls). (23) 

where F'^ is given in (17). 

Corollary 7 Under Nagakami-m fading, a uniformly randomly chosen connected node x G $ has mean 

6(6 + m) 

= U ' (24) 

■ms[d + 1) 

which is 1 + 6/m times the value without fading. 
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Proof: A random connected node x is distributed according to 

/.(a;) = ^ . (25) 
^ ^ EAT 

Without fading, the distribution is s^5x^~^, < x < resulting in an expectation of 5/ {s{6 + 1)). ■ 
For Rayleigh fading, for example, the density is a gamma density with mean 5/s, so the average 
connected node is 1 + (5 times further away than without fading. 

B. Connectivity with retransmissions 

Assuming a block fading network and n transmissions of the same packet, what is the process of nodes 
that receive the packet at least once? 

Corollary 8 In a network with iid block fading, the density of the process of nodes A" that receive at 
least one of n transmissions is 

A"(x) = (1 - F{sxY)cd5x^-^ . (26) 

Proof: This is a straightforward generalization of Prop. 6(a). ■ 
So, in a standard network, the number of connected nodes with n transmissions 

EAT" = I A"(ar)da; = — (*(n + 1) + 7) , (27) 

Jo s 

where ^' is the digamma function (the logarithmic derivative of the gamma function), which grows with 
logn. Alternatively if the threshold Sk for the /c-th transmission is chosen as Sk = si/k, k e [n], the 
expected number of nodes reached increases Unearly with the number of transmissions. 

IV. Broadcasting 

A. Broadcasting reliability 

Proposition 9 For S = 1 and Nakagami-m fading, m G N, the probability that a randomly chosen node 

X G [0, a) can be reached is 



k\ 

k=0 



, (28) 



where s = as. Pm is increasing in m for all s > and converges uniformly to 

lim = min{l, s~^} . (29) 
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Proof: Pm{s) is given by 

= F{sx))dx = C £K!!!£^dx . (30) 



/o Jo r(m) 

For m G N, this is 

— (msx) 
k=Q 

which, after some manipulations, yields 

-1 k 



Pm{s) = 2_^ I exp(-msx) — - — dx , (31) 



Pm{s) = ^ f 1 - — exp(-ms) V V 



ms) 



(32) 



1 - exp(-ms) 2^ jtT^ 

A;=0 



(33) 



The polynomial Pm-\ is the Taylor expansion of order m of (1 — s) exp(ms) at s = (the coefficient 
for s™ is zero). So exp(— m.s)Pm_i(s) = 1 — ,s + 0(5"*+^) from which the limit 1 for s < 1 follows. 
For s> 1, the exponential dominates the polynomial so that their product tends to zero and 1/s remains 
as the Umit. ■ 
The convergence to min{l, s~^} is the expected behavior, since without fading a node is connected if 
it is positioned within [0, 1/s] (s < 1) and for a randomly chosen node in [0, a] for a > 1/s or s > 1, this 
has probabiUty 1/as. So with increasing m, derivatives of higher and higher order become at s = 0. 
From the previous discussion we know that Pm{s) = l + 0(s"^). Calculating the coefficient for s"* yields 

The m-th order Taylor expansion at s = is a lower bound. Upper bounds are obtained by truncating 
the polynomial; a natural choice is the first-order version 1 + (m — l)s to obtain 

^~ r(m + 2) ^'") <i^"»(^) ^min|l,^(l-exp(-ms)(l + (m-l)s))| . (35) 
Using the lower bound, we can establish the following Corollary. 

Corollary 10 (e-reachability.) If 

»»<(n=i±^i^. (36) 

m 

at least a fraction 1 — e of the nodes G [0, a) are connected. In the standard network (specializing to 
m = 1), the sufficient condition is 

as <2e, (37) 
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This follows directly from the lower bound in (35). 
Remarks: 

- For m — >^ oo, the bound (36) is not tight since the RHS converges to 1/e for all positive e (by 
Stirhng's approximation), while the exact condition is as < 1/(1 — e). 

- The sufficient condition (37) is tight (within 7%) for e < 0.1. With pi{as) = (1 - e~"*)/as, the 
condition pi{as) > 1 — e can be solved exactly using the Lambert W function: 

as < >V(-ge~«) + q , where q = . (38) 

A linear approximation yields the same bound as before, while a quadratic expansion yields the 
sufficient condition as < 2e + 4/3e^ which is within 3.9% for e < 0.25. 

B. Broadcast transport sum-distance and capacity 

Assuming the origin a transmits, the set of nodes that receive the message is {xj}. We shall determine 
the broadcast transport sum-distance D, i.e., the expected sum over the all the distances x^^"^ from the 
origin: 

L> = E [ ^ x^/" j (39) 
VxgI- / 

Proposition 11 The broadcast transport sum-distance for Nakagami-m fading is 

5 1 r(m + A) 

Dm = Cd-r . ^ ' ("^0^ 

A (msj" i [m) 

and the (broadcast) fading gain Dm/Doo is 

A^^ 1 r> + A)^ 

Doo r(m) ^ ^ 

Proof: From Campbell's theorem 

/•oo 

= Cd6 / x^/"+^-\l - F{sx))dx , 
Jo 

which equals (40) for Nakagami-m fading. 



(41) 
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Without fading, a node Xj is connected if Xj < 1/s, therefore 

Doo= x^/"A(x)da; (42) 
Jo 

=«1'"=^^4t'"- '''' 

So the fading gain D^/D^ is the A-th moment of f as given in (41). ■ 
Remarks: 

1) The fading gain is independent of the threshold s. Dm oc for all m. It strongly resembles the 
connectivity gain (Prop. 6), the only difference being the parameter A instead of 5. In particular, D„j 

is independent of rn if A = 1. See Remark 3 to Prop. 6 and Fig. 4 for a discussion and visualization 
of the behavior of the gain as a function of m, and A. 

2) For Rayleigh fading (m = 1), Di = Cd6s~^, and the fading gain is r(l + A). For d = a = 2, 

n - 

3) The formula for the broadcast transport sum-distance reminds of an interference expression. Indeed, 
by simply replacing x^/" by x~^, a well-known result on the mean interference is reproduced: 
Assuming each node transmits at unit power, the total interference at the origin is 



E 





which for (5 < 1 diverges due to the lower bound integration bound (i.e., the one or two closest 
nodes) and for (5 ^ 1 diverges due to the upper bound (i.e., the large number of nodes that are far 
away). 

So far, we have ignored the actual rate of transmission R and just used the threshold s for the sum- 
distance. To get to the single-hop broadcast transport capacity C (in bit-meters/s/Hz), we relate the 
(bandwidth-normahzed) rate of transmission R and the threshold shy R = log2(l + s) and define 

C = max{i? • D{2^ - 1)} = max{log2(l + s)D{s)}. (44) 



Let be the broadcast transport sum-distance for s = 1 (see Prop. 11) such that Dj^ = Dj^s ^ 
Proposition 12 For Nakagami-m fading: 

(a) For A G (0, 1), the broadcast transport capacity is achieved for 

>v(-^)+A- 



February 2, 2008 



DRAFT 



17 



The resulting broadcast transport capacity is tightly (within at most 0.13%) lower bounded by 

log 



Cm > %^(A-^ - A) (e^-^-^ - l)"^ . (46) 



(b) For A = 1, 



Cm = ^, (47) 
log 2 



independent of m, and Rgpt = 0. 
(c) For A > 1, the broadcast transport capacity increases without bounds as 0, independent of the 
transmit power 

Proof: 

(a) Djji oc s~^, so Cm oc R{2^ — 1)~^ which, for A ^ 1, has a maximum at i?opt given in (45). The 
lower bound stems from an approximation of i?opt using W(— exp(— 1/A)/A) ^ —A which holds 
since for A = 1, the two expressions are identical, and the derivative of the Lambert W expression 
is smaller than -1 for A < 1. 

(b) For A = 1, Cm increases as the rate is lowered but remains bounded as R ^ 0. The hmit is 
Cdd/\og2. 

(c) For A > 1, R{2^~1)-^ is decreasing with R, and limR^o i?(2^-l)-^ = limR^o(log 2)-'^R^-^ = 
oo. 

■ 

Remarks: 

- The optima for R, s are independent of the type of fading (parameter m). 

- For A < 1, the optimum s is tightly lower bounded by 

sopt ^ exp(A-^ - A) - 1 . (48) 

This is the expression appearing in the bound (46). 

- (c) is also apparent from the expression D(s)log2(l + s), which, for s ^ 0, is approximately 
DmS^~^/\og2. So, the intuition is that in this regime, the gain from reaching additional nodes 
more than offsets the loss in rate. 

- For A = 1/(2 log 2), Sopt = ^opt = 1 and Cm = Dm- This is, however, not the minimum. The 
capacity is minimum around A 0.85, depending slightly on m. 

Fig. 5 depicts the optimum rate as a function of A, together with the lower bound (A~^ — A)/ log 2, 
and Fig. 6 plots the broadcast transport capacity for Rayleigh fading and no fading for a two-dimensional 
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Optimum rate vs. A 
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Fig. 5. Optimum transmission rates for A £ [0.5, 1.0] The optimum rate is 1 for A = 1/(2 log 2) « 0.72. 



network. The range A G [0.5, 1.0] corresponds to a path loss exponent range a G [3, 6]. It can be seen 
that Nakagami fading is harmful. For small values of A, the capacity for Rayleigh fading is about 10% 
smaller. 



C. Optimum broadcasting {superposition coding) 

Assuming that nodes can decode at a rate corresponding to their SNR, the broadcast transport capacity 
(without fading) is 

C = E ^xi/"log2(l + x-i) (49) 

To avoid problems with the singularity of the path loss law at the origin, we replace the log by 1 for 
X < 1. For X > 1, we use the lower bound log2(l + x~^) > Proceeding as in the proof of Prop. 11, 
we obtain 

C>CdS(^^ + J^^ x^-^dx^ , (50) 
which is significantly larger than in the case with single-rate decoding. For A < 1, 



C> 



A(l - A) 



(51) 
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Broadcast transport capacity vs. A 
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Fig. 6. Broadcast transport capacity for d = 2, A G [0.5, 1.0] and m = 1 and m = oo. For A = 1, the capacity is 
27r/(31og2) w 3.02 irrespective of m. For the no fading case, the minimum occurs at A = 1/(2 log 2), where C = 2n/3. 

For A ^ 1, this lower bound and thus C is unbounded, in agreement with the previous result. The only 
difference is that for A = 1, C diverges whereas C is finite. Note that since log2(H-a;~^) < l/(xlog2) 
for X > 1, the lower bound is within a factor log 2 of the correct value. 

If the actual Shannon capacity were considered for nodes that are very close, C would diverge more 
quickly as A ^ (a ^ oo) since the contribution from the nodes within distance one would be: 



A. Maximum transmission distance 

How far can we expect to transmit, i.e., what is the (average) maximum transmission distance M = 



Let X be a uniformly randomly chosen connected node. The pdf fx is given by (25). The distribution 




(52) 



V. Other Applications 



E(max^g^{xV°})? 
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of the maximum xm of a Poisson number of RVs is given by the Gumbel distribution^ 

Fs,^{x) = exp (-EiV(l - Fs,{x)j . (53) 

So, in principle, M = E(x|^") can be calculated. However, even for the standard network, where 
-Fxm(^) = exp(— ^ exp(— sx)), there does not seem to exist a closed-form expression. If the number of 
connected nodes was fixed to ca/s (instead of being Poisson distributed with this mean), we would have 

^xm{^) = (1 ~ e^^s)Cd/s ^ijj^ mean 

ExM = ^(*(y + l)+7) ■ (54) 

Since ^' is concave, this upperbounds the true mean by Jensen's inequality. Finally, we invoke Jensen 
again by replacing E(x|^") by E(xm)^^" to obtain 

^<(;(*(7 + 0+^))"°- 

Without much harm, "^{x) could be replaced by (the slightly larger) log(.T). Even replacing 'I'(.x+ 1) by 
log(a;) still appears to be an upper bound. The bound is quite tight, see Fig. 7. Also compare with Fig. 1, 
where the most distant node is quite exactly 6 units away (s = 0.1). The factor is the bound in 

the non-fading case, so the Rayleigh fading (diversity) gain for the maximum transmission distance is 
roughly log(l/s)^/" which grows without bounds as s — >^ 0. 

B. Probabilistic progress 

In addition to the maximum transmission distance or the distance-rate product, the product distances 
times probability of success may be of interest. Without considering the actual node positions, one may 
want to maximize the continuous probabilistic progress G{x) = max{x^/"P[/^ > sx]}. For the standard 
network with a = 2, this is maximized at a; = l/2s. If there was no fading, the optimum would be 
X = \/l/s. Of course there is no guarantee that there is a node very close to this optimum location. 

Alternatively, define the (discrete) probabilistic progress when transmitting to node i by 



G,- = E f x//" • W > sx,- 1 X 



]) (56) 



We would like to find Zopt = arg maxj Gj. For the standard network, 

n w M/" ( \\ 4 r(i + l/a) 

G. = E(x, ^M-sx^)) = ^^—-^^^ . (57) 

'Note that the Gumbel cdf is not zero at 0"*". This reflects the fact that the number of connected nodes may be zero, in which 
case the maximum transmission distance would be zero. Accordingly the pdf includes a pulse at 0, the term exp(— EiV)(5(a:). 
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Expected max. transmission distance vs. threshold s 
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Fig. 7. Expected maximum transmission distances for the standard two-dimensional network and s € [0.05, 1.00]. For 
comparison, the curve s~^^^ for the non-fading case is also displayed. 



The maximum of Gi cannot be found directly, but since r(i + l/a)/r(i) is very tightly lower bounded 
by i^/" we have 

J jl/a 

Gi ^ , '' \-^u (58) 
which, assuming a continuous parameter i, is maximized at 

iopt = — i — ,] , s\ ■ (59) 
alog(l + ^) 

Note that the same expression for igpt 

would be obtained if Gi was approximated by the factorization 

G'i = E(x//°)P[^i < 1/s]. For the standard network, E(x//") = ^(;+Y/"\ and P[^i < = (vr/(7r+s))\ 
So G[ differs from Gi only by the factor (1 + s/cdY^°' which is independent of i and quite small for 
typical s. 

Now, the question is how to round iopt to iopt- For large s, iopt = 1- For small s, iopt ~ Cd/{ois) so 

^opt = r— 1 (60) 

as 

is a good choice. It can be verified that this is indeed the optimum. The expected distance to this iopt-th 
node is quite exactly l/(as)^/". So in this non-opportunistic setting when reUability matters, Rayleigh 
fading is harmful; it reduces the range of transmissions by a factor q;~^/°. 
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C. Retransmissions and localization 

Proposition 13 (Retransmissions) Consider a network with block Rayleigh fading. The expected number 
of nodes that receive k out of n transmitted packets EiVj^ is 

^^fe=^^^7r^' ke{0,l,...,n}. (61) 



Proof: Let p{x) = 1 — F{sx). The density of nodes that receive k packets out of n transmissions is 
given by 

XUx) = \{x) - P(x))"-^ . (62) 

Plugging in p{x) = exp(— sx) for Rayleigh fading and integrating (62) yields ^NJ^ = A^(R''"). ■ 
Remarks: 

- Interestingly, (61) is independent of n. So, the mean number of nodes that receive k packets does not 
depend on how often the packet was transmitted. 

- Summing over k € [n] reproduces Cor 8. 

- (61) is valid even for /c = since KNq = oo. 

- For the standard networks, the expression simpUfies to EiVj^ = ^, which, when summed over € [n], 
yields (27). 

Let be the position of a randomly chosen node from the nodes that received k out of n packets. 
From Prop. 13, the pdf (normahzed density) is 

f^„{x) = XUx)-^^y ke[n]. (63) 

For the standard network, we have Ex^ = (ns)~^, Varx^ = (ns)~^, and Exf = ^(*(n + 1) + 7), 
which is again related to (27) (division by the constant density c^). 

The densities of the nodes receiving exactly A; of 6 messages is plotted in Fig. 8 for the standard 

network with a = 2. 

This expression permits the evaluation of the contribution that each additional transmission makes to 
the broadcast transport sum-distance and capacity. 

These results can also be applied in localization. If a node receives k out of n transmissions, Ex^ is 
an obvious estimate for its position, and Var x^ for the uncertainty. Alternatively, if the path loss x can 
be measured, then the corresponding node index i{x) can be determined by the ML estimate 

i{x) = arg maxj (x) , (64) 
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X 

Fig. 8. Densities A|(a:) for the standard network with a = 2 (cd = ir) and s = 1. The maximum of the density for fc = n = 6 
is A6(0) = 7r. The dashed curve is the density of the nodes that receive at least 1 packet. NormaUzed by EN^ these densities 
are the pdfs of . 

with the pdf given in Cor. 2. For the standard networks, for example, the ML decision is i{x) = \cd/x] 
since 

i{x) = i ^ ^^x<-^. (65) 
I I — 1 

This is of course related to the fact Exj = i/cd- 

VI. Concluding Remarks 

We have offered a geometric interpretation of fading in wireless networks which is based on a 
point process model that incorporates both geometry and fading. The framework enables analytical 
investigations of the properties of wireless networks and the impact of fading, leading to closed-form 
results that are obtained in a rather convenient manner. 

For Nakagami-m fading, it turns out that the connectivity fading gain is the (5-th moment of the fading 
distribution, while the fading gain in the broadcast transport sum-distance is its A-th moment. A path 
loss exponent larger than the number of dimensions d (d+1 for broadcasting) leads to a negative impact 
of fading. Interestingly, the broadcast transport capacity turns out to be unbounded if A > 1, i.e., if the 
path loss exponent is smaller than d -I- 1. While this result may be of interest for the design of efficient 
broadcasting protocols, it also raises doubts on the validity of transport capacity as a performance metric. 
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Generally, it can be observed that the parameters S and/or A appear ubiquitously in the expressions. 
So the network behavior critically depends on the ratio of the number of dimensions to the path loss 
exponent. 

Other applications considered include the maximum transmission distance, probabilistic progress, and 
the effect of retransmissions. We are convinced that there are many more that will benefit from the 
theoretical foundations laid in this paper. 
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